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In this work, an Effective Field Theory (EFT) incorporating light SU(3)-flavour and heavy quark spin 
symmetry is used to describe charmed meson- antimeson bound states. At Lowest Order (LO), this means 
that only contact range interactions among the heavy meson and antimeson fields are involved. Besides, the 
isospin violating decays of the A (3872) will be used to constrain the interaction between the D and a D* 
mesons in the isovector channel. Finally, assuming that the X(3915) and Y(4140) resonances are D*D* and 
D*D* molecular states, we can determine the four Low Energy Constants (LECs) of the EFT that appear 
at LO and, therefore, the full spectrum of molecular states with isospin J = 0, i and 1. 
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1. Introduction 

Since the discovery of the X(3872) resonance 
many theoretical approaches have been used trying 
to describe it. Even though it contains a cc pair, it 
does not seem to fit in the charmonium spectrum. 
Because of that disagreement, other more exotic 
proposals have been made. Among these, the inter- 
pretation of the X(3872) as a hadronic molecule is 
the most likely so far. Within this assumption, the 
X(3872) would be a bound state of both a charm 
meson and antimeson. However, this hypothesis 
strongly depends on the J PC quantum numbers of 
the resonance, which have not been experimentally 
established. They are either 1 ++ or 2 h (analysis 
done in Ref. Q), of which only 1 ++ is compatible 
with a low-lying S-wave hadronic molecule. 

Apart from the X(3872), many other experimen- 
tal hidden charm resonances have been observed: 
the XYZ states. Many of these states can be new 
candidates for heavy meson- antimeson molecules. 

These resonances are near their threshold. This 
means that the meson and the antimeson are not 
so close to be sensitive to find the details of the in- 
teraction at short distances. So, both of them are 
preserving their individuality and will not probe 
the specific details of the short range interaction 
responsible of their binding. Hence, a scale sepa- 



ration exists and the interaction between these two 
mesons can be described by different EFT's @,|3]. 
Besides, the presence of a heavy quark in the heavy 
mesons imposes that our EFT should be consistent 
with Heavy Quark Spin Symmetry (HQSS), which 
implies certain constraints to the heavy meson- 
antimeson interaction |5[. On the other hand, the 
light quark content of the heavy mesons {q = u,d,s 
in this work) imposes SU(3) flavour symmetry and, 
because of that, molecular states should be classi- 
fied into SU(3) multiplets. 

As a consequence of the symmetries discussed 
(HQSS and SU(3)) flavour symmetry) only four 
parameters are enough to describe the molecular 
states at leading order in the EFT we use. That is, 
we need four data points to predict the full molec- 
ular spectrum. For this purpose we will assume the 
molecular nature of certain XYZ states such as the 
X(3872), X(3915) and Y(4140) and the fourth as- 
sumption will be derived from the analysis of the 
isospin violating branching ratio of the X(3872) de- 
cays into J/'i/jLO and J/ipp. 

2. The EFT Description at Lowest Order 

In this section, the EFT used in this work to 
describe the heavy meson molecules is briefly pre- 
sented. The EFT description must involve pions 
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and heavy meson/antimeson fields and local inter- 
actions among these degrees of freedom that are 
compatible with the known low energy symmetries, 
most notably HQSS and chiral symmetry. Never- 
theless, according to [f|: pion exchanges are weaker 
than naively expected and only enter as a perturba- 
tion at subleading orders. Similarly, coupled chan- 
nel effects turn out to be also sufficiently small to be 
ignored at LO. Hence, at LO, the EFT consists of 
heavy mesons and antimesons interacting through 
a contact range potential, similar to the one in [5|. 
Then, once we have determined our potential, the 
wave functions and observables will be calculated 
in the standard quantum mechanical fashion. For 
example, we can generate bound states by iterating 
the EFT potential in the Schrodinger / Lippmann- 
Schwingcr equation, as previously done in [6|. Sec 
details in 

3. Isospin Symmetry Violation in the 
X(3872) 

The first step is to solve the bound state equation 
for_the A (3872). We consider that the A (3872) is a 
DD* molecule with quantum numbers J PC = 
where we distinguish between the neutral (D°D*°) 
and charged (D + D*~) components of the wave 
function. That is, there are two channels in the 
bound state equation. And we will regularize the 
LO potential (here, and in the following cases when 
not specified) with a Gaussian regulator function 
with a cut-off A varying in the 0.5 — 1.0 GeV range. 

Moreover, the Belle collaboration reported 
the decays of the A (3872) into the (isoscalar) 
J/ , i'7r + 7r~7r and the (isovector) J/^tt + 7t~' chan- 
nels. The latest measurements yield to the isospin 
violating ratio 0: 



Bx = 



T(A(3872) -> J/^tt+tt-tt ) 
T(A(3872) ->• J/fTT+Ti-) 



= 0.8 ±0.3 



which is difficult to accomodate from the theoretical 
point of view. 

In this work, we will assume that this isospin vi- 
olation is caused by the isospin breaking generated 
by the mass difference of the neutral (D°D *) and 
charged (D+D*-) channels in the A(3872), which 
would not have a definite isospin. In this picture, at 
short DD* distances, the X(3872) would be a linear 
combination of / = and 1=1 components whilst 
the J/^3tt and J/^2tt decays would be described 



by an isospin invariant coupling via an intermediate 
p and uj meson (as suggested in [9(). 

The analysis performed by Hanhart et al. in Q, 
leads to the following branching ratio 



Rx 



0.26 



+0.08 
-0.05 ' 



(2) 



M(X ^ J/*w) 

where Bx is translated into a ratio of the decay am- 
plitudes of the A (3872) instead of its corresponding 
decay widths. In our model, this ratio Rx can be 
rewritten as (details in pj) 



Rx 



g P ^xo - *xi 



9u> f 



xo 



(3) 



XI 



being g v = M V [DD*{I = Q, 1) J/* V] with 
V = p/oj and ^xo y ^ xi an average of the neutral 
and charged D*D wave function components in the 
vicinities of the origin. 

So, with the experimental determination of Rx 
and the binding energy of the A (3872) we could 
determine the contact range potential that binds 
the A(3872) if we knew the g p j 'g u ratio,. But the 
9p/ 9u ratio can be determined from the SU(3) re- 
lation g p — g^ = —\plg^ , and, as the OZI (the 
strange quark pair creation is suppressed) rule im- 
plies 3 P ,<7 W 3> <70 , that is, we are ignoring the 
X(3872) decay into J/ip<j>. Thus, we can approx- 
imately consider g p /gu — 1, so the Rx ratio in our 
model is 



Rx 



xo 



A'l 



(4) 



*X0 + *X1 

which only depends on the wave function compo- 
nents in the vicinities of the origin. 

4. The SU(3) and HQSS Partners of the 
A(3872) 

If we determine the value of the counterterms of 
the LO EFT, we will be able to calculate the loca- 
tion of the molecular partners of the X(3872). 

There arc four unknown LEC's. We fix two of 
them from the location of the A (3872) resonance 
and its isospin breaking branching ratio, as ex- 
plained in the previous section. The remaining 
two require the identification of two partners of the 
A (3872): the A (3915) as a 0++ D* D* molecule 
and the F(4140) as a 0++ D*D* molecule, guided 
by its apparently dominant decay into J/^cf). 

Apart from the Rx errors, there is an extra error 
source which has to be taken into account: the ap- 
proximate nature of HQSS. In this EFT, we expand 
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HH 


E (A = 0.5 GeV) 


E (A = 1 GcV) 


0++ 


DD 


3709 ± 10 


3715ti5 


1++ 
1+- 


D*D 
D*D 


Input 
3815 ± 17 


Input 
382lt$ 


0++ 

1+- 

2++ 


D*D* 
D*D* 
D*D* 


Input 
3955 ± 17 
4013 n 9 


Input 

3958^27 
4013^2 





HH 


E (A = 0.5 GcV) 


E (A = 1 GeV) 


0+ 


D+D- 


3835.8l|| 


3837.7^;| 


1+ 


D S D* 


3949 ± 13 


3957tig 


0+ 
1+ 
2+ 


D*D* 
D*D* 
D* S D* 


4056 ± 22 
4091+^ 


4061^ 
4097l 2 g 



Table 1: Predicted masses (in MeV) of the SU(2) isoscalar HQSS 
partners of the X(3872) resonance for two different values of the Gaus- 
sian cutoff. Errors in the predicted masses are obtained by adding in 
quadratures the uncertainties stemming from the two sources of sys- 
tematic errors discussed at the end in Subsect. [4] tf : see discussion in 

a 



Table 2: Predicted masses (in MeV) of the isospinor 
(/ = |) HQSS partners of the X(3872) resonance, 
for two different values of the Gaussian cutoff. The 
meaning of the quoted errors in the table is the same 
as in Tabled 



the QCD lagrangian into powers of (Aqcti/™q) so 
that 



V LO 



yhO 



(1± ^w ); 



(5) 



Taking an approximate value of m c ~ 1.5 GeV 
for the charm quark mass and Aqcd ~ 200 MeV, 
we should expect a 15% violation of HQSS for the 
LO contact range potentials. Since the two error 
sources are independent, the total error will be com- 
puted by adding in quadratures the partial errors. 

4-1- The SU(2) Isoscalar (I — 0) Partners 

In this sector the hidden strange components 
are ignored. We do not take into account particle 
coupled channel effects as they are subleading un- 
less the mass gap between the neutral and charged 
channels is similar in size to the binding energy in 
the isospin symmetric limit. Thus, in the 1 ++ and 
2 ++ channels, we are using a coupled channel po- 
tential. The analysis of this sector was previously 
done in Ref. Q without including strangeness into 
the analysis (SU(2) light quark flavour symmetry) 
and, as there is almost no difference in the reso- 
nances predicted, we conclude that the effect of the 
isospin violation and the inclusion of the quark s 
must be small. The spectrum of molecular states is 
presented in Table [TJ 

4-2. The Isospinor (I = |J Partners 

In this sector, the C-parity of the molecules is 
not well defined, as they are not bound states of a 
heavy meson and its antimeson. However, the for- 
malism is identical to the one in the previous case, 



except for the 1+ D S D* and DD* molecules. The 
D S D* and DD* thresholds are separated by only 
2 MeV and require a coupled channel treatment. In 
this case, we obtain just one single bound state. 
This resonance and its other isospinor partners are 
shown in Table [2j where we have considered only 
the strangeness one states (it would be the same 
for strangeness S = — 1). 

4-3. The Isovector (1=1) Partners 

The potential in the isovector and isospinor sec- 
tor is the same (except in the 1 ++ and 2 ++ 
molecules owing to the isospin violation). For this 
reason, the spectrum in the isovector sector would 
be similar to the isospinor one. The four molecular 
states obtained are listed in Table [3] The other two 
possible states correspond to the isovector partners 
of the V(3872) and A"(4012) resonances where a 
single bound state is obtained. 

4-4- The Hidden Strange Partners 

In this sector, the resonances must also con- 
tain a ss quark- ant iquark pair. The potential de- 
rived for this latter case is the arithmetic mean of 
the isoscalar and isovector one. So if there was 
a bound state in the isoscalar and isovector sec- 
tor, there would probably exist its hidden strange 
partner. On the other hand, considering that the 
V(3872) and AT(4012) molecules have no isovector 
partners, it is very likely that there will not be 
hidden strange partners either with those quantum 
numbers. Therefore, the four states obtained are 
listed in Table H 
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D* S D S 
D* S D S 
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Table 3: Predicted masses (in MeV) of the SU(2) isovector 
HQSS partners of the X(3872) resonance for two different val- 
ues of the Gaussian cutoff. The meaning of the quoted errors 
in the table is the same as in Tabled 



Table 4: Predicted masses (in MeV) of the hidden strange 
isoscalar HQSS partners of the X(3872) resonance for two dif- 
ferent values of the Gaussian cutoff. The meaning of the quoted 
errors is the same as in Tabled 



5. Conclusions 



As a summary, we have established the location 
of 15 molecular partners of the A(3915), Y(4140) 
and AT (3872) states using HQSS and SU(3) flavour 
symmetry, see Tables [T]|4] These predictions have 
a series of uncertainties, being the most important 
the approximate nature of HQSS. According to the 
estimations of these uncertainties, the spectrum of 
these states should be stable but their location can 
vary up to a few tens of MeV in certain cases. 

However, the family of D^D^ states we the- 
orize depends on the assumptions we made: the 
molecular description of the A(3872), AT(3915) and 
Y(4140) resonances. Whilst there is a consensus on 
the molecular nature of the X(3872) resonance, and 
its existence is well established, the situation for the 
A(3915) and Y(4140) are more dubious. Thus not 
all the predictions are probable alike. Predictions 
derived from the A (3872) are supposed to be more 
solid than those depending on the A(3915), which 
in turn are less hypothetical that the ones obtained 
from the Y(4140). In this sense, as mentioned in 
Ref. @] too, the 2 ++ D*D* isoscalar partner of the 
A (3872) is still the most reliable prediction of the 
present work. So, if any new analysis of the XYZ 
resonances provides us with a better molecular can- 
didate than one of our assumptions, it can easily be 
included in this scheme. Thus, we will be able to re- 
place one of the doubtful molecular states we have 
assumed with this new resonance so the new pre- 
dictions are more robust. On the other way, if even- 
tually any of the predictions we have established is 
detected, it could be a proof to the molecular nature 
of these resonances. 
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